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We study the tunneling between two quantum Hall systems, along a quasi 
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i . one-dimensional interface. A detailed analysis relates microscopic parame- 

ters, characterizing the potential barrier, with the effective field theory model 
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for the tunneling. It is shown that the phenomenon of fermion number frac- 
^ . tionalization is expected to occur, either localized in conveniently modulated 

barriers or in the form of free excitations, once lattice effects are taken into 
account. This opens the experimental possibility of an observation of frac- 
tional charges with internal structure, close to the magnetic length scale. The 
coupling of the system to external gauge fields is performed, leading us to the 
exact quantization of the Hall conductivity at the interface. The field theory 
approach is well supported by a numerical diagonalization of the microscopic 
Hamiltonian. 
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I. INTRODUCTION 



The relevance of the edges in the quantum Hall effect was stressed several years ago in a 
seminal paper by Halperin Jl[ . The basic idea is that the existence of gapless excitations at 
the edges, which follows from the general principle of gauge invariance, provides a mechanism 
for the universal character of the quantum Hall effect, or in other words, the presence of 
plateaus for the Hall conductivity, independently of factors such as the degree of disorder 
and the sample geometry After this initial observation, the subject was left untouched for a 
while until the works by Stone M and Wen |3],^j], which showed many interesting connections 
with the quantum field theory of (l+l)-dimensional models , in particular a relation to Kac- 
Moody algebras J5J. One of the consequences of these works is that they may be useful 
as a promising tool for the study of tunneling between quantum Hall systems |J. Such 
a research program is completely supported by the recent advances in the fabrication of 
microstructures, which may give a real oppurtunity for testing ideas originated from (1+1)- 
dimensional quantum field theory. 

Motivated by this possibility, in this paper we will concentrate on the integer quantum 
Hall effect in the simplest situation, viz, the tunneling along a quasi one- dimensional inter- 
face, taking into account that both samples have their first Landau levels completely filled. 
We will make contact here with the theory of fermion number fractionalization which 
has found some applications in the condensed matter physics of polymers |]T0| , |TTj1 , superfluid 
3 He , and recently suggested to play an analogous role at the interface of quantum Hall 
systems |l3|Jl^| . In fact, as we will show, similarly to the case of polyacetylene [ 10| , |TT , [15 ,16 



lattice effects may induce the presence of fractional charges moving along the interface. It is 
interesting, however, to observe that in the present situation the excitations have a peculiar 
internal structure, composed by two "subparticles" , each one localized around one of the 
edges. 

This paper is organized as follows: in section II we present simple arguments which 
give the effective model of tunneling. This will be useful to show what is to be expected 
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from a more rigorous analysis, to be seen in further sections. In section III we build the 
effective model starting from the microscopic definition of an interface. For the sake of a 
more complete analysis, and to study the conductivity, we also consider the presence of 
additional external gauge fields. In section IV we explore some physical consequences of 
the model defined in the previous section, such as charge trapping in specific barriers or via 
an Aharonov-Bohm effect. The Hall conductivity is found to be quantized, regardless the 
presence of localized states or disorder at the interface. In section V we study lattice effects 
in the tunneling, showing that lattice distortions may occur, associated to a gap in the 
fermion spectrum as well as to fractionally charged excitations at the interface. Distortions 
are represented by a complex bosonic field, interacting with the fermion system. In section 
VI we perform a numerical analysis of the problem of charge trapping in a modulated barrier, 
obtaining a very good agreement between the microscopic system and the effective (1+1)- 
dimensional field theory. In section VII we conclude our discussion and point directions 
for future investigations. Finally, in order to make the paper as self-contained as possible, 
there is an appendix on fermion number fractionalization in field theory, where we show in 
detail the computation of the fermionic current for the model of tunneling, according to the 
method of Goldstone and Wilczek [g]. 



II. HEURISTIC MODEL 

We will show how the general form of the effective model of tunneling may be found 
through phenomenological arguments, relying only on a few basic assumptions This 
will provide us with some motivation before taking the more complicated task of a complete 
microscopic analysis. In this section we assume the absence of two-body interactions and 
eletromagnetic perturbations. 

Our problem is to study what happens when we approximate two planar samples, through 
a common plane, taking into account that both of them have their first Landau levels 
completely filled. As one sample gets closer to the other, there will be some tunneling 
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through an approximately one- dimensional interface. A sufficient condition for the presence 
of tunneling is that we have a certain degree of disorder at the interface and that the distance 
between the samples be of the order of the magnetic length £. 

If we consider one of the samples as the square —L/2 < x < L/2 , —L + A < y < A, 
with A ~ £, its edge can be defined by —A < y < A. The system is under the influence of 
a magnetic field B = Bz and its edge may be physically generated via the introduction of 
an electric field E = Ey, which avoids the presence of electrons (at zero temperature) in 
the region y > A. Working in the Landau gauge, A = (—By, 0) , we have eigenf unctions 
localized only in the y direction. According to Stone || we can define the charge density 
operator at the edge as 

/■A 

j(x) = <p + {x,y) <p{x,y) dy , (2.1) 

J -A 

where (p(x, y) is the field operator in second quantization, constructed as a sum of only first 
Landau level states (which is a good approximation to the case of strong magnetic fields). 
We expect the low energy excitations at the edge to be associated to deformations of the 
quantum Hall droplet, contained in —A < y < A [p,P7|. The Fourier expansion of j(x), as 
given by ( |2.1| ), may be recovered from 



H = J i/)t(x)(-ivd x ) M*) dx > ( 2 - 2 ) 

a Hamiltonian of "right-going" chiral fermions. In this expression, v = cE/B is the drift 
velocity of electrons in the sample. From (|2.2|) it follows that if the sign of v is changed, 
the electric current is inverted. In this case we would have a system of "left-going" chiral 
fermions, given by the field %Pl{x). This opposite situation is exactly what occurs at the 
edge of the other sample. Therefore, the transference of electrons between the samples will 
be given by the introduction, in the Hamiltonian, of the tunneling operator tip^tpL (and 
its hermitian conjugate), where t is the amplitude for the tunneling of electrons across the 
interface. We claim, thus, that 



H 



^{-ivd x ) ipp> + i/j£(ivd x ) ip L + tip^ip L + r^i^l dx , (2.3) 



is the Hamiltonian which describes the tunneling process. Taking t = — (ti — U2) 

1^ 



and using the chiral representation for the 7 matrices, 



(2.4) 



7° = ^i 



7 1 = -«72 (2.5) 



7 5 = 7°7 1 = ^3 



where a±, a 2 and are the Pauli matrices, it is a simple matter to show that we can obtain 



( p.3|) from the following Lagrangian 

C = ^(i^do + iv^dx)^ + $(ti + ihi h )^ . (2.6) 

The fact that t may depend on space and time variables, opens interesting experimental 
possibilities, related to the phenomenon of fermion number fractionalization. We will come 
back to this point later. 

III. MICROSCOPIC DERIVATION 

In our phenomenological description, given by the Lagrangian Q2.6|) , the basic input is the 
tunneling amplitude t. A deeper question, therefore, is to ask for a microscopic derivation of 



(2.6), in which the starting point of analysis is the exact Hamiltonian of the system. Only 
in this way we would know how to obtain the tunneling amplitude from a specific potential 
barrier, characterizing the relevant class of microscopic structures for the observation of 
interesting phenomena. 

A. COMPUTATION IN THE ABSENCE OF EXTERNAL GAUGE FIELDS 



Let us consider a system of independent electrons in the two-dimensional (x, y) plane, 
under the influence of a magnetic field B = Bz and confined only in the x direction by the 
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strip | a; | < L/2. We will be interested in the limit L — > oo. A simple model of interface is 
described by the Hamiltonian 

H = ^-[(P x + eBy) 2 + P;]+eV(x,y) . (3.1) 

The first term in the Hamiltonian represents an non-interacting two-dimensional electron 
gas. The second term can be decomposed as the sum of two contributions: 

eV{x,y) = eV 1 {y) + eV 2 {x,y) , (3.2) 

where 

eV l (y) = -Ay 2 ; A > (3.3) 

b 2 



eV 2 {x,y) = g(x)exp f -— I . (3.4) 



The potential V\(y) is a parabolic barrier, which means, as we will see, that for a certain 
range of the chemical potential \x there will be a region \y\ < yo completely free of electrons. 
By varying fi we can have yo G (0, oo). The other piece of V(x), the term V2(x,y), breaks 
translation invariance in the x direction, generating a modulated tunneling amplitude along 
the interface. It is clear that the Hamiltonian ( |3.1p is unbounded from below, but this does 
not present any problem in our approach: we could regularize the potential V\(y) to be 
well-behaved for \y\ ^> £, which would not change the physics of tunneling at the interface. 

When g(x) = 0, we have an exactly soluble model. In this case the time independent 
Schrodinger equation is 

{i [{P* + eBy) 2 + P 2 }- V} <p(x,y) = E<p(x,y) . (3.5) 

Imposing periodic boundary conditions in the x direction, we can look for a solution of the 
form ip(x,y) = exp(ik n x) £(y), where k n = ^p, and n is an integer. The equation for £($/) 
is 

{i [{kn + eBy)2 + p y ] ~ Ay2 } ^ {y) = my) ' (3 ' 6) 



that is 



p2 

y 

2m 



+ 



2m 



A 



2m 



A 



1 (eB) 2 
4 m 2 



(eBf 
2m 



A 



k 2 + — k 2 
K + 2m K 



eB , . 
2^' ; 



£(y) = Ei{y) 



(3.7) 



Therefore, equation (|3.7| ) represents an one-dimensional harmonic oscillator, with eigenfunc- 
tions given by 

1 / a 4 , \ TT / a" 



Cn, P (y) = exp 



2a 2 



H P y- 



(3.8) 



where 



a: 



2m • 



2m 



.4 



2m • 



/ 1 



V2m£ 4 



.4 



(3.9) 



and H p is the Hermite polinomial of order p. The energies are 



E 



n,p 



l\ 1=1 



1 - I ?)" 



(3.10) 



where uj c — is the ciclotron frequency. 

If the external magnetic field B is high enough we can limit the Hilbert space to the 
first Landau level. Therefore, from now on we will restrict the state space to the set of 
normalized wavefunctions 

1/4 ( 



1 



na 2 L 2 



1 



(3.11) 



The electric current associated with a function (p n (x,y) can be immediately computed: 



(jy) = -(P y )=0 

(J») = -(P*-^y) = - 

m i z m 



1-1? 



(3.12) 
(3.13) 



If we give a look at expression ( |3.11| ), we easily recognize an interesting relation between 
the k n space and the real (x, y) space. The wavefunction ip n represents a delocalized state 
in the x direction and a localized state in the y direction. These states are centered in the y 



direction around y n = jrk n . From (|3.10| ) and (|3.11|) , one can see that for each value of the 



energy, there are two states, one centered at +y n and the other at —y n , so that the lesser 
the energy, the greater \y n \. 

We see from Q3.13j ) that n > and n < define, respectively, electrons moving to the 



right and to the left directions along the x axis. We consider a > £, which is equivalent 
of saying that the wavefunctions are spread in the y direction whithin a region of order £. 
This means that two wavefunctions, one with n > and the other with n < will have 
some overlap only if their orbitals are distant from the line y = by a lenght of order £. It 
is, thus, enough to consider b = £ in the expression for V^x,?/), if one wants to produce a 
modulated tunneling amplitude at the interface. 

For each value of the chemical potential \i in the first Landau level, there is a coordinate 
l/o in such a way that all the states contained in the region \y\ < yo will be empty ones. The 
value of yo is a function of \i and enters in our model as a phenomenological parameter. Since 
a > £, it is adequate to take y ~ £ so that there is some overlap between wavefunctions 
situated at the opposite sides of the interface. 



Acording to ( |3.11| ), the state centered at y (or — y ) defines a value of k n (or k_ n = —k n ) 
given by k = yo£ 2 /a 4 . This state has the quantum number n = kL/2n. All we need to do 
from now on is to find a theory for the modes near k and —k. 

Using the wavefunctions we found for the case of g(x) = 0, given by ( p.ll| ), we may write 
the second quantized field operator as 

= a n Vn{x, y) + J2 a n <Pn(x, y) . (3.14) 

n<0 n>0 

In this representation the Hamiltonian becomes 

H e ff = J J H(p(x,y) dydx , (3.15) 

where H is the first quantized Hamiltonian Q3.ip . 

In order to identify the filled states of our model with a "Dirac sea" in the effective 
theory, it is necessary to make some redefinitions of the operators and a^. Let us make 
the following transformation 
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,R 



,R 



(3.16) 
(3.17) 



The equation ( 3.14 ) now becomes 



(3.18) 



n<n 



n>—n 



In the "large box" limit, the set of modes k n becomes dense and we can define the continuum 
theory through 

L 



E 



2tt\ 1 / 2 



L 



a k 



We, thus, obtain 



7r) [ - a k Vk-k dk + (^-) f _ a L k Lp kVk dk . 

In J Jk<k \Z7T/ Jk>-k 



■ dk 

Z7T 

ifk = ¥n(L) 1/2 



1 \ 2 



(3.19) 
(3.20) 



(3.21) 



In order to completely identify the filled states with the "Dirac sea", it is necessary to 
redefine the energy too, setting to zero the energy of the modes k and —k. In this way, we 
shift ( |3~T0D to 



Eh 



(k 2 - k 2 ) 



2m 



1 - 



a 



(k 2 



2k 



(3.22) 



where 



m 



a 



> . 



(3.23) 



Substituting (|3~2T|) into (gl5|) , we get 



H. 



eff 



k<k 
1 



dkE k _ k a k a k + 



k>-k 



dkE k+k a k a k + 



+ 



— / <s dkdk' dx |exp i(—k + k' + 2k)x g(x)c\(k, k') a k + a k , + H.c.j + 

l r + 

— / _ dkdk! dxexp [i(—k + k')x] g(x)c2(k : k') a k a k , + 
2n Jk,k'<k 

dkdk' dxexp[i(—k + k')x] g(x)c2{— k, — k') a k + a k , , (3.24) 



2-7T Jk,k'>-k 
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where 
ct(k, k') 



7r cr 



and 

c 2 (k, k') 



dyexp 



dyexp 



2a 2 



2a 2 



y-^(k-k)) + (y-SLtf + k) 



The above gaussian integrals can be computed exactly. We have 



ci(k, k') 



b 2 



or 



+ b 2 



exp 



a 



2(a 2 + b 2 ) 



( fc + fc / ) 2 -((A;-A;) 2 + (A; / + fc) 2 ) 



y_ 

b 2 



(3.25) 



yl 
i) 2 



(3.26) 



(3.27) 



c 2 (M0 



b 2 



exp 



cr 
2F 



2(a 2 + 6 2 ) 



(k + k' - 2fc) 2 - ((jfc - fc) 2 + (k' - k) 2 ) 



(3.28) 



As we defined before, b ~ a 
ci(k,k') 
c 2 {k, k') 



so that we can write 







f£ 2 ■ 


(£ 


exp< 


[2 [ 






f£ 2 ■ 




exp< 


[2 I 
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(k + k?) 2 -((k-k) 2 + (V + ky) 



~(Jfe + k' - 2k) 2 - ((jfc - kf + (k' - k) 2 ) 



(3.29) 
(3.30) 



The existence of tunneling depends essentially on the function g(x). If, for example 
g(x) = constant, the term in ( p. 24 ) representing the tunneling between electrons R and L 
vanishes. This occurs because in momentum space g(k) = 5(k), whereas in the tunneling 
process each electron changes its momentum by approximately 2k or — 2k. Therefore, an 
interesting class of functions is given by g{x) = e~ 2lkx f(x) + c.c, where f(x) is dominated 
in momentum space by modes uj « k. That is to say, writing f(x) = J exp(iu>x)f(uj) du, 
we are considering functions centered at a; = 0, with support in the interval Au << k. 
Substituting the proposed form for g(x) into (|3.24j), we obtain 
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'1, +k a k a k + 



H ef f= I _dkE k _ k afa k H + [ _dkE k 

Jk<k Jk>-k 

dkduj[f(u) Ci(k, k — u) a k a k _ u + H.c.] + 



+ 



u>—k<k<k 



- I dkdco[f(u;) c 2 (k, k + 2k - u) af a^ +2S _^ + H.cl + 

I k<k,k<-k+u) 



+ / dkduj[f(Lu) c 2 (-k, -k-2k + uj) a k + a L 

I k>~k,k>~3k+u 



k+2k-u 



+ H.c. 



where, according to ( |3.29| ) and ( |3.30| ) 

c\(k, k — uj) 



2 j exp r 2 



I W V 



(3.31) 



(3.32) 



C2(/c, + 2/c — uj) — cx(k,k — uj) 



(3.33) 



C2{~k, —k — 2k + uj) 



1\3 



2j eXP "2 



*-2 +2 *J + 2 (*~2 



(3.34) 



The expression (|3.32|) shows that for k ~ ~ we have Ci(/c, fc — cj) ~ _i= e — . The 
Hamiltonian ( |3.31| ) describes the complete system, and we want to study only the modes 
associated with tunneling, which are given approximately by \k\ < k. We have to select in 
( p.31| ) only those degrees of freedom involved in the tunneling process, retaining the most 
relevant couplings. Therefore, we obtain from (|3.31|) , a new effective Hamiltonian 



H. 



lfr i _dkE k _ k a k + a k + I _dkE k+k a k + a k + (-) exp[-(H 2 )]x 

[fc[<fc J\k\<k \Z/ 



x / dkduj[f(uj) of ajL* + H.c] + c 2 (-2k, 0) f _ dkdu[f{u) a^ 2 - k+ X + H.c.] + 

J\k\<k J\k\<k 

(3.35) 



+ c 2 (0,-2£;) Jkdu[f{u) a L k + a L k+2k 



i\k\<h 
H.c] + H' 



\k\<k 



where 



H' = / _ _dkE k _ k a k + a k l + _ _dkE k+k a k + a 

J\k+2k\<k J\k-2k\<k 



(3.36) 



We can retain in H\fu only the first three terms. In fact, using ( |3.33j ) and ( |3.34| ) we have 



c 2 (-2k,0) = c 2 (0,-2k) 



V=2 



-3(ke) 2 



(3.37) 
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That is, the ratio between C2(— 2k, 0) and the factor -^e ^ 2 in the tunneling amplitude is 
e -2(fcf) 2 _ ]\j OW; jf ^ ^ we have e~ 2 ^ 2 ~ e -2 ~ 10 _1 , which indeed shows that we can 



keep only the first three terms of ( |3.35| ), to investigate the tunneling. 

The expressions for E k _ k and E k+k , can be linearized around k — 0. Using Q3.22| ) we get 



E k _ k = kv 



E, 



k+k 



—kv 



(3.38) 



Taking into account these approximations, we obtain the effective Hamiltonian 



H eff 



/ dk kv a k + a k — / _dk kv a k + a k + 

J\k\<k J\k\<k 



\k\<k J\k\<k 

+ c f _ dkduj[f{uj) af"a£_ w + H.c] 

J\k\<k 



where 



< |-V exp[-Ck£f 



1\5 



2 J ex P 



.'70 



Defining now 



i) L {x) 



2vr 
1 



dkexp(ikx) a k 
dkexp(ikx) a k , 



we see that H^jf is, in coordinate space, 
dx 



H eff 



d d 

^R(-iv—)ilfR + $i( iv fa)i>L + cf(x) ipRipL + cf*(x) xplipR 



which agrees with ( |2.3| ). 



(3.39) 



(3.40) 



(3.41) 
(3.42) 



(3.43) 



B. INTRODUCTION OF GAUGE FIELDS 

We will obtain now the effective Hamiltonian for the interface, taking into account the 
presence of an external gauge field a^. We have to consider the more general microscopic 
Hamiltonian 
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H 



2m 



(P x + eB y - ea{) 2 + (P y - ea 2 ) 2 + eV(x, y) + ea . 



(3.44) 



We can write 



H = U + H n U 



(3.45) 



where 



H ° = 2m~ \-( Px + eBy)2 + ^ + tV[X ' ^ + 6a ° 



(3.46) 



and U is the unitary operator 



U = exp — ie / a ■ dx' 



(3.47) 



In this expression c represents a path in the (x, y) plane, given by 



c : 



y — , < x' < x 
x — , <y' <y 



(3.48) 



Let us supose that a„ is a small static field with slow variations in the magnetic length 



scale. From H , equation Q3.46Q , we are led, according to our previous computations, to the 
effective Hamiltonian 



u 1 

U 0,eff 



dx 







+ dx ea (x, y = 0) feV'i? + ^1 V> 



(3.49) 



The second quantized unitary operator U takes the form 



U ~ 1 — ie / + (x, y) 



ai(x',0) dx' + / a 2 (x,y') dy' 



(f>(x, y) dxdy . (3.50) 



Using ( |3.21| ) and neglecting terms similar to a§ + a%,, we obtain 



U = 1 



ie 
2^ 



/ _ dkdk' dx exp[i(—k + k') x] / aAx', 0) dx' c 2 (k, k')[af + ay + 

Jk,k'<k ' Jo 

dkdk' dx exp[i(—k + k') x] Cz(x, k, k') af + a§ + 



+ *\*W - IT . , , , 
Z7T Jk,k'<k 



ie 

2n Jk,k'<-k 



dkdk' dx exp[i(—k + k') x] ai(x, k, k') a k + a^ 



(3.51) 
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where 



c 3 (x, k, k') 



dy / a 2 {x,y') dy' exp 



2a 2 



a 



v-«r(*-*) + 



a 



(3.52) 



Ci{x, k, k') 



( 1 



V7ra 2 



(/// / a 2 (x,y') dy'exp J 



af 



(* + *)) + 



a 



(3.53) 



and £2 (A;, A;') = lim^oo C2(fc, k'), with c 2 (fc, A;') given by ( 13.281) . We must retain in equation 
(P.51J) only the modes associated with tunneling, \k\ < k. Approximating a 2 (x, y) by a 2 (x, 0) 
and substituting k = k' = in the expressions for c 2 , c 3 and c 4 , we get 



U J = 1 



2ll J\k\,\k'\<k 



dkdk'dxexp[i(—k + k')x] / ai(x',0) dx'laf a5 + 

'|<fc Jo 



L+ Li * e 

27T 7|fc|,|fc'|<fc 

+ Z/oafc +a fc') a 2 (z,0) 



dkdk'dx exp[i(—k + k')x] ■ (—y a§ + a§ + 



(3.54) 



or alternatively, using the definitions ( p.41| - |3~4~2j ), 

U 1 = 1 — ie J dx ipRipR (—yoa 2 (x, 0) + J ai(x', 0) dxj + 
+ i^ti'L (yo^ix, 0) + / ai(x',0) ete' 



(3.55) 



To compute H*** = U 1 Hq^jj-U 1 , it is important to know the following operator products: 



i) U 1 Vfl(») ^ = 
iP~r{x) + ie dx' 



iI)r{x')iI) R (x') -y a 2 (x', 0) + / ai(x", 0)dx'' Vis fa) 



■0a + ^ e ( ~~ yo a 2(%, 0) + y Oifa', 0) dx' ) ^sfa) 



exp 



( — yoa 2 (x, 0) + y ciifa', 0) da/ 



V^fa) 



(3.56) 



and, in the same way, 
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ii) U I+ ipti x ) U 1 = exp ie (yoa 2 (x,0) + J a\(x',0) dx' 



tpi(x) . 



(3.57) 



Therefore, we get, using the above relations, 



U I+ 

= ^ + R 
-41 











dx 



dx 



U 1 



d d f f x 

dx ~ iC ~dx V ~ y ° a2 ( x > °) + J 01 ( x ' °) dx ' 

d 



d of f x 1 1 

dx ~ ie dx \ y ° a2 ( x > °) + J ai ^ x '' °) dx 

dx 



fa 



(3.58) 



where D x = ^ — ieax(x, 0), and also 



iv) U 1 i/JriPlU 1 = exp(-2iea 2 (x, 0) y ) ip^tpL 



(3.59) 



From (|3.58| ) and (|3.59 ) we obtain 



Hi f} = U I+ H^ ff U T = J dxi-i^vD^ - c^ifr + if 2l ">) exp(2tea 2 y 0l ^ + 



The gauge invariant Lagrangian associated with this Hamiltonian is 



(3.60) 



C = #(7° A) + v^-Dx) ^ + c^{f x + ^ 27 5 ) exp(2*ea 22 /o7 5 )^ + eva 2 y d^ . (3.61) 

It is not dificult to understand why the above Lagrangian is the correct generalization 
of ( |2.6|) , which takes into account the presence of gauge fields. The first attempt, in order 
to implement gauge invariance, would be to replace, in ( |2.6|) , 7 M <9 M (we are now considering 
v — 1) by the covariant derivative, Jf) = 7°(<9o + iea ) + 7 1 (<9i + iea{). This is, however, 
only a partial answer. It is necessary to consider that the operators i/)r(x) and iPl{x) 
create hole states in different positions of the two-dimensional plane. In this way, if a gauge 
transformation is performed, ipR and ipi will be multiplied by different phase factors. Besides 
that, one additional requirement is that the effective model of tunneling be gauge invariant 
at a classical level. This condition is obtained from the decoupling betweeen the modes at 
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the interface and the bulk degrees of freedom. Considering the gauge transformation (below 
A* = 0, 1) 



a fi a fi + d^a and a 2 — > a 2 + <9 2 a| y= o 



(3.62) 



associated to 



^ and ^ — > ^>ea;p 



i —7 + ea 



(3.63) 



where = e9 2 «l2/=o2yo is the phase factor implied by the physical distinction between ip R 
and ipL, as fields defined in edges separated by 2|/ , a simple check shows us that (|3.61| ) in 
fact satisfies all the above conditions. 

We have not considered, in the microscopic derivation, the presence of two-body inter- 
actions. Indeed, as far as the wavefunctions at the edges are delocalized in the x direction, 
their low electric charge densities make the Coulomb repulsion generate logarithmic effects 
on the spectrum [IB]. The Coulomb potential may be relevant, however, in situations where 
we have localized states, induced by specific configurations of the tunneling amplitude. 



IV. CHARGE TRAPPING AND CONDUCTIVITY AT THE INTERFACE 



Let us explore some of the consequences of the effective model of tunneling, given by 
( |3.61| ). We will consider first the case do = ai = a 2 = 0. This model may exhibit the 
phenomenon of charge fractionalization 0, intrinsically related to the global behaviour of the 
external fields fi(x) and / 2 (x). A gradient expansion of the fermion current, computed 
assuming fi and / 2 as slowly variating fields in space-time, yields || (see appendix) 

1 



_^ /u/ ^ab fa 9y fb 



tan " I 



(4.1) 



(27T) |/| 2 (27T) 

For a configuration given by fi = e, e — ► and / 2 = v4tanh(Ax), we obtain, from the 
expression for a total charge +1/2 or —1/2, localized near x — 0, when e — > + or 

e — > 0~, respectively. The role of /1 is limited to providing a regularization of ( O ). The two 
possibilities for the total charge are associated with a zero mode in the spectrum, which by 
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its turn implies a doubly degenerate state: the occupied zero mode, with charge 1/2 and the 
unoccupied one, with charge —1/2. According to the discussion of section III, this profile of 
/ corresponds to the modulating potential 

g(x) = 2Atanh(Ax) sin(2A;x) . (4.2) 

This potential could, in principle, be manufactured with the techniques used in the fabrica- 
tion of microstructures. A difficulty is that the spatial period of (|4.2|) occurs at not yet very 
manageable scales (~ lOOipA), by present day technology. However, it is possible that in 
disordered interfaces (which must contain, in Fourier space, modes close to 2k), modulated 
in adequate scales by a function qualitatively similar to tanh(Az), we have a more pratical 
way to observe the occurrence of fermion number fractionalization. 

Let us now investigate the introduction of perturbing gauge fields in the system. There is 
a curious Aharonov-Bohm effect, which shows how a magnetic flux may control the amount 
of charge localized at the interface. In order to see it, we consider a = 0, / = const, and a 
magnetic flux $, confined in the interior of a very thin imaginary solenoid crossing the plane 
at the point (x, y) = (0,0). In this way, the exponential factor which appears in ( |3.61| ) is 
exp(— ze$7 5 ), for x = 0~ and exp(ie$7 5 ), for x = + . According to (f4.1|), this discontinuity 
will induce an accumulation of charge $/7r at x = 0. 

It is important to remark that equation ( |4.1|) , while giving a good evaluation of the total 
amount of charge trapped by variations of the tunneling amplitude, must be replaced by a 
sharper computation if we want to know the size of the region where most of the charge is 
confined. This may be found through exact solutions of the Dirac equation, in special cases. 
In a "soliton" profile, as / = iAtanh(Xx), the zero mode may be explicitely computed [19| 
and the localization region inferred to be 

( 2xt> \ 

5 ~ lim rx ... ., , . 4.3 

x ^°° Vic U\f{v)\dy) 

We also note that in the situations we have been discussing, the charge concetrated at 
the interface will have an interesting internal structure: the relation between Fourier space 

17 



and coordinate space (the y direction) will make any charge distribution have two equivalent 
and disjoint pieces, localized at the edges R and L. The numerical analysis of section VI will 
show this effect very clearly. 

Since we know how the system is coupled to gauge fields, we can compute the conductivity 
at the interface, as a response to small electric fields. 

The current density which crosses the interface, j±, flowing from one edge to the other, 
may be calculated through 

% 51nZ[a 2 ] i 5\n.Z[a 2 + a] 



where a = a(x,y) and 



2y 5a 2 2y 



Z = Di^Dipexp 



5a 



i I d 2 xC 



(4.4) 



a=0 



eff 



(4.5) 



is the generating functional of model ( |3.61| ). We can now use the chiral anomaly PU] of 



model ( |3.61| ) to find 



Z[a 2 + a] = Z[a 2 ] ■ exp 



In 



(4.6) 



where F^ u = d^a? — d u a 11 . In this way, from ( |4.6| ) and (^4]), we obtain 



3± = ~2^ Ex 



(4.7) 



Above, E x is the x component of the electric field. We may also study the influence of a 
small electric field, pointing in the y direction. In this case, we have 



a = a 1 = and a 2 = E v t 



(4.8) 



where E y represents the y component of the electric field. Applying ( flip to this problem, 

we get 



(J 1 ), = ^-e 10 d (-2ea 2 y ) = ^-E y 2y 

Z7T Z7T 



(4.9) 



Therefore, the above relations tell us that the response to external electric fields is given (in 
the usual units) by the following conductivity tensor 



a 



a 



XI) 



\ 



-a 







with <j xv = — , 
y h 



(4.10) 



showing the quantization of the Hall conductivity at the interface. It is interesting to note 
that this result is strongly related to the imposition of gauge invariance and is independent on 
the specific configuration of the tunneling amplitude, even when it is associated to localized 
states. The above argument may be regarded as an exact version, worked out for a particular 
problem, of the fundamental works of Halperin [I]] and Laughlin [ 2Tj] , , since here we do not 
have to take any averages over magnetic fluxes. 



V. LATTICE EFFECTS 

The quantum Hall effect is observed in approximately two-dimensional electron systems, 
confined at the interface of semiconductor devices as Si-SiC>2 or GaAs-Gai_ x Al x As. The 
latter has been preferred in recent years due to its high-mobility parameters, allowing for 



very precise measurements of conductivity and other physical quantities |22 . 

Let us consider for a study of lattice effects, the two-dimensional electron gas as having 
a certain thickness 2Lj_, experimentally found to be around 50(pA, and interacting with a 
three-dimensional atomic lattice of macroscopic size, in fact much larger than any relevant 
characteristic length in the process of tunneling. The lattice structure is relatively complex. 
The GaAs lattice, for instance, has a zincblende structure, consisting of two interpenetrating 
fee sublattices, one of Ga and the other of As atoms. The inter-atomic distance, s , is close 
to 5(fA. We will study the effects of the lattice on the tunneling, through an approach 



analogous to that of Takayama et al. [|15j for the case of polyacetylene. 

The fortunate fact that the magnetic length is nearly lOOtpA in the quantum Hall effect, 
means that we can, in order to estimate couplings, simplify the discussion assuming the 
lattice to have a simple cubic structure composed of interlinked atoms through an harmonic 
potential. Therefore, the lattice potential may be expressed as 
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v (x) = [ V a(x - Xi) - ^(x^djVoix - Xi)\ , (5.1) 

i 

a sum of single atomic potentials centered at sites Xi = s{n\x + n 2 y + n s z), with nj, n 2 , n 3 
integers, where we consider distortions as represented by £(afi), in a linear approximation. 
We will be interested in static distortions of the lattice, so that we can write the total 
Hamiltonian of the system (quasi one-dimensional interface + lattice) as 

H = Hi ff + -K (Ite) - fe)) 2 + / d 3 x^(x)<P(x)V(x) , (5.2) 

4 <M> J 

where #<f// is given by ( |3.43| ), K ~ 10~ 6 (v9t4)~ 3 is the force constant associated 
to distortions, and <f>(x) is defined from ( |3.21| ), according to (f>(x,y) — > 4>(x,y,z) = 
(2/(7rL^_)) 4 exp (— (z/L±) 2 ^ (f)(x,y) , which corresponds to three-dimensional wavefunctions 
confined in \z\ < Lj_. The sum in (|5.2j ) is carried over nearest neighbors. We are also 



assuming the absence of additional external gauge fields. 

A convenient approximation is to replace sums over sites in ( |5.1|) and (|5.2| ) by integrals. 



We find, then 



V(x) = i J d 3 xV (x) + dPSid&ixJVoix - Xi) (5.3) 

and 

H = H{ ff + — J tfZd&ffid&ix) + J d*x4> + {x)(t>{x)V{x) . (5.4) 

As an estimate, we can write Vq(x) = —47rZe 2 s 2 5 3 (x), where e 2 = 1/137 is the fine structure 
constant and Z ~ 2 is an effective atomic number. Substituting the above quantities in 
V(x), equation ( [5.3D , we obtain 

AitZp 2 Ait 

V (x) = _ —Ze 2 d^(x) . (5.5) 

The first term on the RHS of ( |5.5|) is a constant and may be shifted to zero. As we noticed 
in section III, only certain Fourier components of V(x) will be relevant in the tunneling. 
This amounts to considering a new set of fields, a(x) and ip(x), given by 
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£(x) ~ a(x) cos(2fcx + ip(x)) 



(5.6) 



From (|5"H| ) we get, neglecting gradients of a and ip, 

dj£j(x) = —2k<r 1 (x) sin(2fcr + f{x)) + cos(2fcr + (p(x))dj(Tj(x) + 
—(jj(x) sin(2A;x + <p(x))djip(x) — —2kai(x) sin(2fcc + f{x)) ■ 



(5.7) 



Therefore, we see that only a\ must be considered in our analysis. On the other side, the 



quadratic term in £ in the Hamiltonian ( |5.4| ) turns out to be (restricted to a\ 



K 



— / d 3 xdj (&i cos(2A;x + ip)\ dj (&i cos(2/cx + <p) 



K 



d A x 



(2k) 2 af sin^(2fcr + tp) + (Vaif cos 2 (2fcx + <p) + 



+ (Vf) 2 sin 2 (2fcx + ip) — dj(Ti(2k5ij + djip) cos(2fc + ip) sin(2fcx + ip) 
K 



4s 



(Fx 



(2k) 2 o{ + (VcnY + (VcpY 



(5.8) 



where, assuming that o\ and <p are dominated by slow modes, we have substituted the above 
trigonometric functions by their averaged values. 

Defining the complex field a = 0\ exp(if) and using the approximations ( |5.7| ) and (p.8|), 
we find, for the total Hamiltonian of the system, 



H = Hi n + ^-lSx 



(2k) 2 \a\ 2 + (d ] a) + (d J a) 



+ 



+ 



4nZe 2 ki 



d 3 x<p + (x)(j)(x) exp(2ikx)a(x) + H.c. 



(5.9) 



It is useful to work in Fourier space, through 

1 



a(x) 



d qexp(iq • x)a(q) 



(5.10) 



(2vr) 3 / 2 

We now define a variational ansatz for the Hamiltonian fl5.9|) : let us consider a class of fields 
&(q), parametrized by Ai, A 2 > 0, according to 



5(g x ,0,0) if \qy\<& x ll , |g,|<A 2 /Lj 
0, otherwise. 



(5.11) 
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This means , roughly, that we are taking the distortion field cr(x) as diferent from zero 
only in a certain neighborhood of the interface, given, in Fourier space by the two varia- 



tional parameters Ai and A2. Substituting (|5.11| ) in ( |5.9| ) we find, after straightforward 
computations, the Hamiltonian 



where 



H = H* + / ,/,- 



eff s£L l 



(d x a(x)) + (d x cr(x)) +ti\<t(x)\' 



+ 



+ 



ic/(Ai,A 2 ) J dxiPt(x)ip R (x)a(x)+R.c. 



(5.12) 



(2tt) 1 /2 



dqexp(iqx)a(q, 0, 0) 



(5.13) 



A 2 A 2 



(5.14) 



and 



/(Ai,A 5 



/ £2 \ a 2/l± 

f% exp — -gj / dg 2 exp 

M/e \ 4 J J~a 2 /l ± 



(5.15) 



(5.16) 



In order to study the variational problem, we will consider the path-integral formalism, 
defining the generating functional 



Z = J DipDipDa+Daexp (iS[ip, a}) 



(5.17) 



where 



S[ip, a] = J (fx^itpi^do + v^d^ip + ip(t 1 + it 2 j 5 - cf{A 1 , A 2 )\a\ exp(-^ 7 5 )) ip + 
KA 1 A 2 



slL\ 



(d l( j) + (<V)+r/M 



(5.18) 



We will simplify our discussion , assuming that ti,t 2 = 0, or in other words, that we have 
a very clean interface, without any degree of disorder or modulating potentials. We can 
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integrate over the fermion fields in (|5.17|) , neglecting variations of |cr| and ip. This may be 
considered as the first term in a gradient expansion of the fermion determinant. We get, 
then, 



J eff 



Da + Da exp (iS e ff[a\) 



(5.19) 



with 



S, 



eff 



(5.20) 



where 



Veff(W\ 



(c/(A 1 ,A 2 )| t r|)- 
2ttv 



In 



' 2cf(A 1 ,A 2 )\a\ 
vk 



+ 1 



(5.21) 



In the computation of (p.21[) it is important to consider the presence of a cutoff at k in 



the fermion theory. The variational strategy is to find extremes of (|5.20|) , in the space of 
configurations of <j(x) and also in the space of parameters Ai and A 2 . Let us perform this 
analysis in two steps: first, we consider a fixed pair (A 1; A 2 ) in order to find the (constant) 
field a(Ai, A 2 ) which extremizes S e ff. Second, we look for extremes of S e ff [a (A 1; A 2 )], in 
the space (A 1; A 2 ). 

We obtain, in the first step, the gap in the fermion spectrum, 

TlvKA^T] ( 1 \ 



\cf (Ai, A 2 ) a\ = vkexp 



-1 - 



c 2 s£L ± \f(A u A 2 ) 
and the effective action, evaluated for a (A 1; A 2 ), 



(5.22) 



Seff [^(Ai,A 2 ) 



KAiA 2 r] 



2s£L A 

EA t A 2 (FA\ + GA 2 2 + H) 
(/(Ai,A 2 )) 2 



\<j\ 



exp 



where 



P _ Kv*k 2 F _ i r 



JAjAg {FA\ + gA| + if) 
(/(Ai,A 2 )) 2 



1 H = {2k)\l 2?:rK 



(5.23) 



(5.24) 



8s£L ± c 2 ' * U 2 ' " 3L\ ' " ' * c 2 s£L_ 

It is readily seen that the second step in the variational analysis, dS e ff /dA\ = dS e ff/dA 2 
0, leads to 
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<9Ai, 



A X A 2 {FA 2 + GA 2 + H) 



(5.25) 



TAiAg {FA\ + GA 2 + H) _ 



(5.26) 



(/(Ai,A 2 )) 

or 

, 2 + G£ 2 
(/(Ai,A 2 )) 

In order to study the above equations, we have to define k and a (t> depends on a; see 
equation ( |3.2ij| )). Let us, then, consider k = 3/(2£), corresponding to an interface with the 
edges separated by ~ 3£, and a/£ = 1.02. With these definitions, we get v ~ 5 x 10" 5 , 
c ~ ((27r) 1/2 /6) x 10~ 4 (ipA)~ 2 and I~4x 10~ 5 . A quick inspection shows that both (|5.25|) 



and (|5.26|) have solutions and that the minimum of S e ff[a(Ai, A 2 )] is obtained from ( ^.26 ). 



This equation has, in fact, many solutions, which determine a curve in the (A 1; A 2 ) plane. 
An estimate yields an isotropic solution A x ~ A 2 ~ 10. If we look at ( 5.11|) , we see that the 



degenerecence in the solutions of ( |5.26| ) could mean a "torsion" in the lattice displacements, 
if the parameters Ai and A 2 were non-trivially dependent on q x . We assume, however, that 
the isotropic solution represents a mean field for the fluctuations of Ai and A 2 . Anyway, the 
gap, as given by ( |5.22| ) does not depend on the degenerecence of Ai and A 2 : substituting 
( p.26| ) into ( |5.22| ) , we find two possible "vacua" , 

cf (A U A 2 ) a ~± lO- 7 ^)- 1 , (5.27) 

physically measurable as a gap at the interface. 

We can, in the same way, find non-trivial solutions of the Euler-Lagrange equations for 
the complex field a. We will have, in general, solitons which interpolate the vacuum values 
of < a >. The transition region is given by the square root of the ratio between the kinectic 
and mass term coefficients in the action ( |5.20| ). This quantity is r^ 1 ! 2 . Using Aj, A 2 ~ 10, 



we have rf 1 ! 2 ~ 10y?A Since a changes its sign in this specific configuration, we will have 
solitons carrying charges it 1/2 (compared to the electron charge) propagating along the 
interface. Using equation ( |4.3| ) we find that the charge will be spread, along the interface, in 
a region of lenght ~ 500(/>A Here, as in the case of polyacetylene, we can also have soliton- 
antisoliton (polarons) excitations. One can suppose that the soliton states would be found 
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as midgap excitations, but since the gap is estimated to be small, it would be hard to observe 
any kind of soliton production threshold through variations of the sample temperature. We 
point, however, that there is a mechanism, related to Coulomb repulsion, which raises the 
soliton energy up to more clearly observable scales. The argument is as follows. As we have 
already noticed, the edges are associated to different wave- numbers of the fields ipR an d i>L- 
This means that all the charge- density profiles of these excitations will be symmetrically 
displaced at opposite sides of the interface. In this way, the Coulomb repulsion will raise 
the soliton energy by ~ (e/4) 2 /(2£) ~ 3 x 10~ 6 ((pA)~ 1 , which is close to the gap between 
Landau levels (~ 4 x 10~ 6 (<^y4)~ 1 ). We see, therefore, that the Coulomb interaction, which 
was playing a minor role in the theory so far, has an important participation in the soliton 
spectrum. 

Another interesting point is the possibility of an enhanced electric conductivity at the 
interface, via soliton excitations. The same phenomenon was conjectured to be present in the 
polyacetylene, but in view of the relatively small polymer filaments, the solitons probably 
do not contribute directly to the conductivity ^,^4j]. In our case, however, the interface's 
length may be constructed many orders of magnitude larger than the magnetic length, in 
such way that we could hope the solitons to be relevant in the conductivity process. 



VI. NUMERICAL ANALISYS 

In order to test the accuracy of the approximations made in sec. Ill we performed a 
numerical investigation of the model of tunneling defined by relation (3.1). 

Let us consider the matrix elements of (3.1) in the basis of functions ip n (x,y) , as given 
by (3.11), and put it in the folowing form 

< n\H\m >=< n\Ho\m > + < n\Hj\m > , (6-1) 

where 

Ho = ^[(P x + eBy) 2 + P*] + eV 1 {y) (6.2) 
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H I = eV 2 {x,y). (6.3) 

The elementary calculation of < n|/f | m > leads to a diagonal matrix with its elements 
given by (3.10), with p = 0. In order to evaluate < n\Hi\m > we must choose a specific 
modulating potential barrier. Since we want to explore the occurrence of fermion number 
fractionalization, we may take a function f(x) with the same asymptotic behaviour as the 
modulating potential considered in sec IV or that of the predicted solitons of sec V. Some 
numerical improvement is obtained from the consideration of 

f{x) = - \e-ig sgn(x) (l - e- A|x| )l . (6.4) 



In ( |6.4| ) we have considered A = (5£) _1 and g = 4/m£ 2 . We have performed the numerical 
diagonalization of ( |6.1| ) considering L/i ~ 900, k = 2/£ and a "smooth" interface, deter- 
mined by [aji — 1) ~ 10~ 4 . The eigenvalue distribution is shown in fig.l. In this figure 
we show the gap due to the presence of Vi{x, y). The order of magnitude of the gap agrees 
with the value obtained from ( |3.43| ), A ~ eg. The two points inside the gap correspond 
to the "zero mode" and to a spurious wavefunction (associated to the periodic boundary 
conditions) . 

In order to obtain the charge distribution along the interface we must evaluate 

\M*,y)\ 2 = Y^WYcftpfayYipfay) , (6.5) 



where is the i component of the n th eigenvector of ( |6.1| ). The result for the localized 
state (zero mode) is shown in fig. 2a. On the other hand, a typical delocalized state in the 
band, far away from the gap, is depicted In Fig. 2b. Note that in the localized state we 
have peaks at both sides of the interface. The integral of the two peaks is close to 1/2, so 
that at each side of the interface there is an accumulated charge of 1/4, measured in units 
of the electron charge. In the numerical computation the total charge is not exactly 1/2, in 
view of finite size effects, as was numerically verified. 

It is interesting to analyse the "vacuum" structure of the theory. For this aim it is enough 
to integrate expression ( |6.5|) in the y direction. We define, then, a projected charge density, 
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|?/; n (x)| 2 = J^dy \i/j n (x, y)\ 2 . In the field theory context, the vacuum charge density is 
evaluated by means of 



calculation. The numerical result, shown in fig 3, agrees with the field theory expression. In 
this computation one observes that the states in the upper band cancel the states in the lower 
band, in such a way that the charge density of the "vacuum" turns out to be determined 
only by the zero mode. This fact suggests that the approximation of linearizing the energy 
near the Fermi energy is indeed a very good aproximation. This result also confirms that our 
matrix Hamiltonian is adequate to obtain the physics at the interface, without interferences 
from the bulk. Therefore, the equivalence of the two calculations presented above is a clear 
evidence that the "Dirac sea" of the field theory model reproduces accurately the completely 
filled lower band. 



We studied the tunneling across a quasi one-dimensional interface in the integer quantum 
Hall effect. The particular form of the electron spectrum at the edges allowed a mapping 
from the microscopic definition of the interface to a relativistic (1 + l)-dimensional quantum 
field theory. Gauge fields and lattice effects were considered in the description. Regarding 
the coupling to gauge fields, the Hall conductivity was found to be quantized, independently 
of the possible induction of localized states by a non-uniform tunneling amplitude. We also 
obtained a peculiar Aharonov-Bohm effect, which shows the influence of magnetic fluxes on 
the charge concetrated at the interface. The study of interactions between edge excitations 
and a three-dimensional lattice showed us a natural mechanism for the generation of fraction- 
ally charged solitons propagating along the interface. They are associated to topologically 
stable solutions of the Euler-Lagrange equations for a complex scalar field. We point that 




(6.6) 




lower band but also the upper band is considered in the 



VII. CONCLUSION 
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these excitations may contribute strongly to the a xx component of the conductivity tensor. 
A numerical test supported the field theory approximations in the case of charge trapping 
in a modulated barrier. 

We believe that an experimental investigation of the above predicted phenomena is cru- 
cial for a further development of the theory, in the sense of a more accurate quantitative 
description. Anyway, there are a certain number of extensions of the present work which 
may motivate future studies. These are related to different microscopic definitions of the 
interface, a consideration of higher Landau levels, the introduction of spin in the formalism 
(which may be relevant to the case of magnetic fields of lower intensity) and a derivation of 
lattice effects, taking into account exact crystal structures (a Monte Carlo numerical analysis 
with pseudopotentials could provide some information on lattice distortions and solitons). 
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APPENDIX A: FRACTIONAL FERMION NUMBER 



In order to make this paper self-contained, we will show in this appendix how the fermion 



current (|4.1|) may be obtained. We will compute it through the adiabatic method, as outlined 
by Goldstone and Wilczek ||. A more rigorous approach may be found in ref. f9j. 
Let us consider the following fermionic Lagrangian in (1 + 1) dimensions 

C = ^ + cij(f 1 + il 5 f 2 )ij , (Al) 

where fx and f% are classical external fields, and the 7 matrices are defined in ( |2.5| ). This 
model may exhibit the phenomenon of fermion number fractionalization, according to the 
topology of the fields fi and f 2 . 
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The model is invariant under rotations in the (fi, f 2 ) plane. This is related to global 
chiral transformations. That is, considering 



VK/i + ^7 5 / 2 ) i> = |/|^exp 



., -l / f2\ 5 

ztan - 7 



(A2) 



where \f\ = {fl + f 2 )^, we see that a rotation in the plane of coordinates / 2 ) by an 
angle gives tan _1 (/ 2 //i) tan _1 (/ 2 //i) + @, which can be absorved by a global chiral 
transformation 



if) -> exp ( -- 67 5 ) V 



■0 — >■ ?/;exp ( — - 07 5 



(A3) 
(A4) 



The existence of this symmetry will be important to establish the topological nature of (J 1 '). 

Let xo be a space-time point such that fi(xo) 7^ and f 2 (xo) = 0. If we consider fi 
and f 2 as slowly variating fields, we can calculate (J M (x )) using a "mass" —cfi(x ) and an 
interaction term £j = icip^ 5 f 2 ip , as can be seen from ( |A1| ). We may write 



(J»(x )) = (0mx ) 7 ^(xo)|0) = lim(O|T0(x o ) 7^(^o + e)|0) 

e o< 



(A5) 



That is , defining x\ = xq and x 2 = xq + e, 



(J /i (xo)) = Tr[A^(x 1 ,x 2 )] 



(A6) 



where A^,Jxi,x 2 ) = 7^ o .(O|T-0 (T (xi) ipp(x 2 )\0) and we are suppressing the limit simbols to 
simplify the notation. 

Expanding ( |A6| ) in a perturbative series, we find, up to the first order in c 



Tr 



A M (xi,x 2 ) = A#(xi,x 2 ) + cjd 2 x A$(xi,x) f 2 (x) 7 5 7^A(^(x,x 2 ) 



(A7) 



where Ag(xi,;r 2 ) = 7 M So(a;i, x 2 ) = 7 M (0|T^(a;i) i>(x2)\0)\f 2 =o- Now, since 



A$(xi,x) =7 m (2tt) 2 y c/ 2 giexp[2gi(xi - x)] S (qi) 
f2(%) = / d 2 kexp(ikx) f 2 (k) , 



(A8) 
(A9) 
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where So(q) — (jq — m) 1 , we obtain 



(J»(x )) = Tr 



C7 A 



/" d 2 q 1 d 2 kd 2 q 2 d 2 x S (?i) / 2 (A;) 7 5 S te) exp[igi(xi -a;)] x 



(2^) 4 

x exp(i/ca;) exp[2g 2 (:r ~ ^2)] 
Integrating in x and qi and substituting g 2 by q, we find 



(A10) 



(J" (a*)) = Tr 



07^ 
x exp(ikxi) 



J d 2 qd 2 k S (k + q) 7 s / 2 (£;) S'o(g) exp[?g(:ri - x 2 )] 



x 



(All) 



In the limit Xi — > x 2 = x , we have 



( J"(x )) = Tr 



07^ 
(2^ 



d qd k S (k + q) 7 / 2 (fc) S'o(g) exp(ikc ) 



(A12) 



Substituting the expansion Sq(A; + g) = <So(<z) + k l, -^So(q) in the above expression, we get 



(J"(x )> = Tr 



+ Tr 



+ 



^ / A 5 (?) 7 5 5 (?) ^0) 

erf 2 {xo) J d 2 q JL So (q) 7 5 <? (g) 

(jq + m)(2q u ) 



(2vr) 2 
icd"f 2 (x ) 

(2^) 2 
icd u f 2 (x ) 



Tr |y y d 2 

Tr[ 7 /i 7 I/ 7 5 ] y d 2 g 



dg 1 

7* 



(g 2 — m 2 ) (g 2 — m?) 



(g 2 — m 2 ) 



m 



(2tt) 2 ^ ' L ' ' ' J J " * (g 2 - m 2 ) 2 ' 
Since Tr[^^^ u rf] = 2e fiu , we obtain 



(A13) 



/A 



(2tt) 2 J ^ (q 2 - m 2 ) 2 
Performing a Wick rotation , q — > iq , the above integral yields 



(A14) 



% \ d 2 q 



m 



(q 2 + m?y 



2%i / dq 



mq 



(g 2 + m 2 y 

■ r a d ( 
—mm / dq 



<ig y g 2 + m? 



m 



(A15) 



In this way, using m = —cfi(x ), we find 
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{■n*o)> 



i 



d„f 2 (x ) 



(A16) 



(2rr) /i(x q ) 

As we mentioned, the theory has chiral invariance. This means that we must write 
the current (J M (a;o)) in a chiral invariant way. In other words, it must be invariant under 
rotations in the (fi, / 2 ) plane. Therefore, we are led to 

-iff* 



I JHfrrW — 1 e P» fa 9 U fb _ 1 nv r, 

{J [X)) ~ (27T) P ~ (27T) " 



tan 



(A17) 



It is interesting to note that from a perturbative calculation and taking into account the 
chiral symmetry of the model, it was possible to find the non-perturbative expression (|A17|) . 

The dependence of (J M ) with the topology of the fields fi and f 2 may be clearly seen by 
considering, for example, the total charge 



Q 



(J°(x))dx 



-A tan- 1 (fl 



1 

2^ 



°° -1 / fi 

di tan — 

-oo \/l 



dx 



(A18) 
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FIGURES 



FIG. 1. 

The energy as a function of the quantum number k. The unit of energy is 10~ 4 ((pA)~ 1 
while k is measured in units of 2k. The dots inside the gap represent localized states. They 
were depicted there only to help in the visualization of their energies. 

FIG. 2. 

The charge density at the interface, \if)(x, y)\ 2 , and their level curves. Fig. 2a) shows the 
profile of the charge trapped in the potential barrier. Fig. 2b) shows the charge density of 
a delocalized state. 

FIG. 3. 

The projected charge density at the interface, obtained from the integration of the charge 
density in the y direction. The total charge is close to —1/2 (empty zero mode). The x 
coordinate is measured in units of the magnetic length and all the wavefunctions considered 
in this computation were normalized to 1 in the interval —200 < x < 200. 
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